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ABSTRACT 
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Given an n-th order spline basis it will be proven that the differentiation matrix is 
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periodic. Furthermore, it will be shown that spline-based bases generate a class of 
compact finite difference schemes. 
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1 Introduction 


The use of wavelets as a basis set for the numerical solution of partial differential 
equations (PDE’s) is a very active area of research. Wavelets, Daubechies or Spline- 
based, provide a very convenient structure for dividing information according to the 
frequency of the information at a given location. This splitting of the data by pro- 
jecting onto a wavelet basis seems ideal for the solution of nonlinear PDE’s where 
shocks or high frequencies might exist in a small portion of the domain, whereas in 
the remainder of the domain the solution might be comprised of very low frequency 
information. 

In this paper the goal is to understand the character of a spline-based wavelet 
numerical method. Numerical results for such a method can be found in [12] and 
[6]. This understanding will be found through the differentiation matrix. In [4] 
the differentiation matrix was found for Daubechies wavelets with periodic boundary 
conditions, and a very high degree of differentiation accuracy known as superconver- 
gence was proven to exist. Furthermore, it was seen that a Daubechies-based wavelet 
method corresponds to a finite difference method with grid refinement in regions 
where small scale data is present. In the more complicated scenario where periodic 
boundary conditions are no longer assumed, it was seen in [5] that the superconver- 
gence is lost. 

In this paper the differentiation matrix for a spline-based wavelet basis with peri- 
odic and nonperiodic boundary conditions will be constructed. It will be proven that 
a spline-based wavelet basis, also, displays superconvergence when periodic boundary 
conditions are assumed. Furthermore, the superconvergence is lost when periodicity 
is no longer assumed. In addition, it will be seen that spline-based wavelet methods 
generate a class of compact finite difference schemes. 

Recall that the differentiation matrix is comprised of three matrices, V = C DC : 

• The first matrix C is the quadrature matrix which maps from point values of a 


1 



function to the approximate scaling function coefficients at the finest scale. 

• The second matrix D maps from the scaling function coefficients of a function 
to the scaling function coefficients of the derivative of the same function. For 
convenience, I will henceforth refer to this matrix as the ‘derivative projection’ 
matrix. 

• In this paper the third matrix is the inverse of the first matrix: C = C~ l . 

This paper is organized as follows: 

1. Introduction 

2. Definitions and Framework: 

Aliasing is reviewed. 

3. Multiple Scale Wavelet Decompositions: 

Wavelet methods use information at many scales, but all properties of wavelet- 
based differentiation are dictated by differentiation in the finest scale approxi- 
mation subspace, Vo. 

4. Differentiation in a Wavelet Subspace: 

Here we begin to develop the relations needed to build the differentiation matrix. 

5. Theoretical Accuracy: 

In this section the theoretical accuracy will be proved. For a spline-based 
wavelet basis using an n-th order spline and assuming periodic boundary con- 
ditions, differentiation is accurate to order 2n + 2. 

6. Equations for Finite Dimensions: 

In finite dimensions all wavelet operations can be written in matrix form. In 
this section, the matrix equations will be derived which will hold for arbitrary 
boundary conditions. 
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7. Examples with Periodic Boundary Conditions: 

Explicit differentiation matrices for the n = 1 and n = 3 spline-based wavelets 
will be given for periodic boundary conditions. 

8. Examples with Nonperiodic Boundary Conditions: 

Again, explicit examples are given. The most important point to note in this 
section is the loss of superconvergence at the boundaries. 

9. Loss of Superconvergence: 

The very important issue of the loss of superconvergence when nonperiodic 
boundary conditions are imposed is explored. 

10. Conclusion: 

The role of aliasing is noted in determining the order of accuracy of V. It was 
shown in [4] that Daubechies-based wavelet methods have a strong relationship 
with ordinary finite difference methods. In this paper, a strong relationship 
between spline-based wavelet methods and compact finite difference schemes has 
been exposed. Furthermore, superconvergence is lost when boundary conditions 
other than periodic boundary conditions are imposed. 
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2 Definitions and Framework 


Please note that the framework, definitions, and even much of the notation in this 
paper come from a paper by M. Unser and A. Aldroubi, see [10]- 

2.1 Fourier Transforms 

Let us set once and for all the definitions of the Fourier transform and the semi- 
discrete Fourier transform used in this paper. The Fourier transform of f(x ) £ L 2 (R) 
is, 

A fOO 

f(0 = / f(x)e' ix dx, ( 1 ) 

J — OO 

and its inverse is, 

1 r°° a 

f(x) = Lj (i)e di - (2) 

The semi-discrete Fourier transform (SDFT) of /*, £ 1 2 {Z) is, 

OO 

/(f) = E e**/*, (3) 

k— — oo 

and its inverse is, 

h = t ( 4 ) 

2.2 B-Splines 

All B-splines in this paper are central B-splines: centered at 0. The B-spline of 
order n is denoted by /3 n (x). All splines of order n, denoted by s n (x), are a linear 
combination of B-splines of order n, see [7] and [8]: 

s"(*) = £c*/r(x-fc), ( 5 ) 

k 

and the space spanned by splines of order n is denoted by S n . The B-spline fl°(x) is 
the box function, which is 1 for x £ [—1/2, 1/2] and zero otherwise. All B-splines of 
higher order are generated from fl°(x): 

p n (x) = p° * 0 n ~\x), ( 6 ) 
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where denotes the convolution operator. Other useful forms of this equation are, 


= (7) 

and 

/ oo 

/? n (y)/3 ?l (* - y) d v, ( 8 ) 

-co 

and finally, 

(! 2n +'( x )=<p\f3 n (x--)>, (9) 

where for /, g £ L 2 (R ), 

/ oo 

f(x)g(x)dx. (10) 

-OO 

We will, also, need the samples of B-splines at the integers: 

b n {k) = 0 n (x)\„ k . (11) 

In order to use B-splines to construct a multiresolution analysis it is necessary 
that contractions and expansions of n-tli order B-splines also be contained in S n . 
This is only true when n is odd, which is, therefore, assumed for this paper. 

2*3 Fourier Transforms of B-Splines 

From equation (6) we see that the Fourier transforms of B-splines are generated 
recursively. That is, the Fourier transform of a convolution is the product of the 
Fourier transforms of the convolved functions. Therefore, the Fourier transform of a 
B-spline is, 

r(0 = «»c n+1 (0, (12) 

where sinc({) = sin i^ and is the Fourier transform of 

2.4 Aliasing and Fourier Transforms 

The relationship between the Fourier transform and the semi-discrete Fourier trans- 
form and the role of aliasing is at the heart of a proof that will be given later in this 
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paper. Therefore, I will reconstruct the straightforward algebra leading to aliasing 
using notation defined above. 

Let be the samples of f3 n (x) as defined above: 

1 f°° * 

K = /?"(*) = — j_" P*(t)e-* k d£. (13) 

Break the integral up as follows: 

1 00 /•(2r+l)jr A 

K = E / / nOe~' ik dt, (14) 

27r r^co^r-lfir 

which is equivalent to, 

l ^ 

K = ir- E / /3 n K+2ir)e- i '‘^, (15) 

Z7r r=-oo 

where I have used the fact that 1 = e~ t27rr/c . If we compare this expression with 
definition (4) we see that, 

oo 

H0 = E /?"(£ + 27T7-), (16) 

r=— oo 

which illustrates the aliasing of frequencies £ + 2 Trr onto frequency £. 

2.5 Spline- Based Wavelet Bases 

As mentioned above, the subspace of n-th order splines, S n , is spanned by n-th order 
B-splines. The following theorem will define the scaling functions for a spline-based 
wavelet basis, see [11]: 

Theorem: The set of functions {<f> n (x — k) : k E Z} with 

OO 

#-(*)= e rhinx-k), in) 

k=—oo 

is a basis of S n provided that the sequence {p} is an invertible convolution operator 
from / 2 into itself. 

The sequence { p } is chosen such that the scaling function satisfies whatever prop- 
erties are specified. That is, one can specify that the scaling functions and wavelets 
be orthogonal which places restrictions on {p}. There are many ways to define the 
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scaling function, but in this paper it will be seen that in the calculation of the dif- 
ferentiation matrix that the sequence {p} will ‘divide out’ producing a result which 
holds for all spline-based wavelet bases produced using the current framework. 

Throughout the paper I will use the notation <j> k (x) to denote translation of <j> n (x) 
by k: 

r k (x) = <f> n (x-k). (18) 

Define the dual of <f > n , <j> n , as 

OO 

r(x) = Y r k p n (x-k), (19) 

k — — oo 

where the sequence {r} is chosen such that 

S k =<M k >. ( 20 ) 

Later the relationship between {r} and { p } will be made precise. However, as men- 
tioned above, {p} and {r} will have no affect on the differentiation matrix. 

The wavelet is defined as, 

oo 

il> n {x)= Y <lkP n {x-k ), (21) 

k— — oo 

and its dual is 

oo 

r(x)= y s k nx-k), ( 22 ) 

k=—oo 

where ft is the ‘B-spline wavelet 1 . For this paper this B-spline wavelet is not needed 
and will not be dealt with, see [10]. Furthermore, as with the scaling function, the 
sequences {< 7 } and {s} defining the wavelet and its dual will be of no consequence in 
the calculation of the differentiation matrix. 

This completes the definitions and framework. 
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3 Multiple Scale Wavelet Decompositions 


One of the underlying strengths of a wavelet decomposition is that it is possible to 
easily decompose a function, or signal, into its components at various scales. To 
illustrate, suppose that we are restricted to a finite number of dimensions, say d. 
Denote the finest scale approximation subspace by Vo, and the coefficients of the 
expansion in Vo by s. A usual wavelet decomposition would appear as, 


V 0 = W l eW 2 (B ... © VVj, © Vj, (23) 

where, as is customary, Wj denotes the wavelet subspace at scale j. A function 
projected onto the above subspaces would, therefore, be represented by d coefficients, 
say c. Under the assumption of finite dimensionality, there will be an invertible matrix 
mapping from s to c: 

c= Ps. 


That is, if g(x ) € Vo with the expansion, 


g(x) = 

i=0 


and one projects into Vj and Wj by, 


d/ 2-1 d - 1 

Pv 1 g(x) = 53 13 

j=0 i=0 


and 


(24) 


d/2-1 rf-l 

rw,j(*) = E E^? < 

j = 0 1=0 

then the matrix form of the above linear transformations is denoted P : 


(25) 


(26) 


' ^ =<#,#>, (27) 

and 

Pd/2+j,i =< >» (28) 

for i = 0, ...,d — 1 and j = 0, ..., d/2 — 1. 


8 



Suppose, now, that we have found the matrix D mapping from s to i, where s 
represents the coefficients of the approximation in Vo of the derivative of a function. 
Of course, no one would choose to work with wavelets and stop the decomposition 
in the subspace Vo- The matrix D does, however, characterize differentiation in any 
combination of wavelet subspaces, and the coefficients c can be found from, 

l=PDP~ x c. (29) 

To restate, D dictates the character of wavelet differentiation. The primary con- 
cern in this paper is to understand D. Therefore, the remainder of this paper will 
not mention wavelet decompositions at multiple scales. 
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4 Theoretical Accuracy 


The accuracy of the spline based differentiation matrix can be proved by observing 
the behavior of the semi-discrete Fourier transform of 

d q =< <^>o? <j>q > (30) 

near £ = 0. It will be shown that the Taylor series of d(£) about £ = 0 is of the form, 

d(() = if + 0« 2 “«) + (31) 

and this implies accuracy of order 2 n + 2. 

4.1 The Semi-Discrete Fourier Transform of {d} 

As above, define d q as, 

dq =< <f>0i(f>q > , ( 32 ) 

and recall from the previous section that this is, 

dq = 5 2p‘ r m J P n ( x - m)d n (x - l — q)dx. (33) 

/,m 

Again, note that as of now all summations and integrals are from — oo to oo allowing 
one to shift summations and integrals without changing the limits, 

dq = ^ ZPi r m J (3 n (x)p n (x + m-l- q)dx. ( 34 ) 

l } m ** 

Let us rename this integral as, 

= J fi n (x)/3 n (x + m — l — q)dx , (35) 

which yields the following expression for d q , 

dq — ) ^ Pf m • (36) 

Break apart the double summation, 

dq — ( 37 ) 

I m 
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and shift the index 1, 


d q = Y1 Pl-1 r m a l-m ■ ( 38 ) 

l m 

Let the superscript “ f ” denote time reversal of a sequence to get, 


dq = ^2p' q -l'^2 r m a l-m- (39) 

l rn 

We see that the right-hand side is a double convolution. We can, therefore, use the 
convolution theorem to get, 

d(() = (40) 


Recall from the previous section that the following expression for r(£) was found: 


r«) = (p'(()i 2n+1 (())-‘- («) 


Using this expression for r(£) in the expression for d(£) we get, 

d(() = x a( - ^ — . (42) 

b 2n+1 (0 

The next two subsections wil derive a(£) and 6(£). After these two expressions 
are found the proof of accuracy will begin. 


4.2 The Semi-Discrete Fourier Transform of {a} 

Define the projection of the derivative of a B-spline onto a B-spline as, 

/ oo , 

nv)r(y - *)iy- («) 

-OO 

Recalling that, 

/«) = •{/«), (44) 

and from the convolution theorem, see the introduction, we get the Fourier transform 
of o-(x), 

a(0 = i(sinc 2n+2 (£). (45) 
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The sequence a p is defined as the samples of a(x) at the integers, 

a P = J f3 n (x)fi n (x - p)dx. (46) 

The semi-discrete Fourier transform of a v is, due to aliasing (see the introduction), 

OO 

a(0 = i ^2 (£ + r)sinc 2n+2 (( + 2irr). (47) 

r — — oo 

4.3 The Semi-Discrete Fourier Transform of {6} 

Recall that the samples of the B-spline f3 2n+ *(x) are found from, 

b 2n+1 = j /3 n (x)(3 n (x - k)dx. (48) 

We know the Fourier transform of /? 2n+1 (x) from the introduction: 

/? 2n+1 (0 = sinc 2n+2 ({). (49) 

Also, from the introduction we know how to get the semi-discrete Fourier transform 
of b 2 k n+l from /? 2n+1 (0, 

OO 

& 2n+1 (0 = sinc 2n+2 (( + 2irq). (50) 

g=— oo 

The semi-discrete Fourier transforms a(£) and b(() will now be used in the proof 
of accuracy. 

4.4 Proof of Accuracy 

The proof of superconvergence begins by first proving two lemmas. 

Lemma A: The series, 

OO 

MO = s inc 2n+2 ({ + 2nq), (51) 

g= — oo 

is an even function of £ and has a Taylor series about £ = 0 of the form, 

■MO = 1 + + C 2 O -f .... (52) 
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Proof: Using the definition of a sine function we get, 


“ sin 2 " +2 K + 2 nq) 

« + 2»,) 2 “+ 2 ' 

We can simplify the numerator since sin(£) is periodic with period 27T, 


*i({) = 


£ 


sin 2n+2 (Q 
(£ + 2Trq) 2n+2 


To show that Si(£) is even first let q = 0 to get, 


( 53 ) 


(54) 


= 0 - !/6f + 1/120T* + ...) 2 ” +2 , 

which is even. Now, sum the terms for q = — 1 and q = 1, 

... = sjn^K) . S 'n 2n+2 (Q 

(f + 2jr) 2 ”+ 2 (( - 2jr) 2n+2 ’ 


to get, 


s({) = s(-0- 


(55) 


(56) 


(57) 


Matching every term with positive q with every term with negative q we see that 
si(£) is an even function. 

We can see from equation (55) that the first few terms in the Taylor series when 
q — 0 are 1 + c\( 2 + c< 2 ^ 4 + .... On the other hand, when q ^ 0 each term in the 
summation is of the form, 

sin 2n+2 (0 
(£ + 2irq) 2n+2 




or 


U + 27r<7) 2 ” +2 P(0 = sin 2n+2 (0- 


Using the series expansions of these expressions we get, 


(59) 


(a, + « 2 { + + ...)P(() =((- ( 3 /6 + « s /120 + ...) 2 " +2 (60) 


or 


(c*i + CV 2 ^ + — £ 2n+2 + h.o.t. (61) 
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This implies that, 


P(0 = 0(( 2n+2 ) + h.o.t. 


(62) 


That is, all terms in Si(0 with q / 0 contribute terms 0(£ 2n+2 ), and therefore, all 
low order terms come from equation (55). This yields the desired result. // 


Lemma B: 

The series, 

OO 

5 2(0 = rsinc 2n+2 (t + 2 irr), (63) 

T — — OO 

is an odd function of £ and has a Taylor series about £ = 0 of the form, 

MO = <M 2n+3 + .... (64) 


Proof: First of all, due to the 2n periodicity of sin(0> MO becomes, 

00 sin 2n+2 f^l 

<»> 

Note that the r = 0 term in this sum makes no contribution to 52(0- Now, note that 
each term in the sum has the following form, 




sin 2n+2 ( 0 

r ({ + 2*r) 2 ’>+ 2 ’ 


( 66 ) 


and, consequently, the function g(£, r) defined as follows is an odd function of £: 


Bit, r) = /(0 r) + /(0 -r) = -(/(-*, r) + /(-£, -r)) (67) 


This implies that S2(0 is an odd function of £ since, 

OO 

MO = 5^0(00- (68) 

r= 1 
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Therefore, the Taylor series of 52(0 about £ = 0 must contain only odd powers of £. 
Consider, now, a typical term of 5 2 (^) for r / 0, 


p«) = 


(£ + 27 rr) 2n + 2 


Rewrite this to get, 


(£ + 27rr) 2n+2 ,P(£) = sin 2n+2 (£). 


Note that the binomial on the left-hand side contains all powers of £, and that 
sin 2n + 2 (£) contains only even powers of £ starting with the term £ 2n+2 : 

(<* + cit + c 2 e + -)p( o = f “ +2 + d 1 e n+> + ... (7i > 

We see now that the first term of P(£) must be ^ 2n+2 /co, but this term will also 
multiply the term in the binomial expression requiring that P(£) have as a second 
term e*i£ 2n+3 where c*i = —ci/(cq). The point is that there exist a nonzero term with 
power 2n + 3 in the polynomial P(£). 

That is, in general the series expansion of the function /(£, r) about £ = 0 contains 
all powers p of £ with p > 2n + 2, and when added to /(£, — r) the even powers of £ 
must cancel leaving as the first term the term with the power 2 n + 3. That is, 


* K ) = a { 2 "+ 3 + ... 


This is the desired result. / / 


Theorem: The semi-discrete Fourier transform of the sequence d * has a series 
expansion about £ = 0 of the following form: 


<*({) = if + 0(f 2 “+ 3 ) + ... 


Proof: Recall that the semi-discrete Fourier transform of {d} is, 


d(0 = T 


b 2n + 2 {() 
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i/ 1 \ _ i E£L-oo(£ + 27rr)smc 2n+2 (£ + 2?rr) 
< ~' V ' E^-oo s *' nc2 " +2 (£ + 2 t rq) 


Now, substitute in the previously derived expressions for a(£) and 6 2n+2 (£) to get, 

(74) 

VS T ^i'H) 

Multiply both sides by the denominator to get, 

oo oo 

d{() sinc 2n+2 {( + 2nq) = i (£ + 27rr)smc 2n+2 (£ + 27rr), (75) 

q=—oo r =— oo 

and combine similar terms to get, 

oo OO 

(d(Q — i() X sinc 2n+2 (£ + 2irq) = 2w ^ rsinc 2n+2 (£ + 27rr). (76) 

Q— oo r =— CO 

Now, apply the results from Lemmas A and B to get, 


(<i(0 - >0(1 + c,{ 2 + ...) = af ! "+ 3 + 


This equation implies the desired result, 


(77) 


J(() =i( + 0 (( ! "«) + 


This completes the proof. // 


(78) 


This proves that the sequence {d} differentiates with accuracy 2n + 2, and, conse- 
quently, the differentiation matrix containing the periodic version of {d} as its rows 
will differentiate with accuracy 2n -f 2. 

The next section will provide explicit examples of the differentiation matrices for 
the first and third order B-splines. Also, the accuracy will be shown computationally. 
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5 Equations for Finite Dimensions 


The equations that will be derived in this section will hold for any spline-based wavelet 
basis defined on a finite-dimensional domain. This includes the equations for both 
periodic and nonperiodic boundary conditions. Note that in this section all summa- 
tions are from 0 to d — 1, and consequently the summation limits will often not be 
shown. 


5.1 Quadrature 

Let g(x ) denote the projection of f(x) 6 L 2 (R) in Vo: 

g(x) = PvJ(x) = (79) 

/=0 

where 

$i =< f,&> - ( 80 ) 

Note that in this section the notation b m (x ) will be used to denote the B-spline 
basis functions. For periodic boundary conditions one could have simply used fS(x — 
m), but for nonperiodic boundary conditions the basis functions near the boundary 
are truncated B-splines. The notation used here is, therefore, general enough for both 
periodic and nonperiodic problems. That is, for a nonperiodic domain bo(x) will be 
only the right-hand portion with support width 1 of the B-spline being used. Now, 
using the B-spline expansion of the scaling function in equation (79) we get, 

g{ x ) = ( 81 ) 

l rn 

Combine the summations, 

g{ x ) ~ ^ > &lPl t 7tibm[ x )i ( 8i ^) 

and sample g(x) at the integers to get, 


S , (k) — ) ^ $ iPl km ( ^* ) ? 

l,m 


( 83 ) 
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or simply, 


Group the terms as, 


9 k — ^ ^ SlPl,mbm,k' 

l,m 


9k — , PLm^m t k\ s h 

l m 

and rename the inner group as 7 to get, 


9 k = 'Yjli'kSi. 

I 

Now, transpose 7 , 

9k = J2 
/ 

and utilize matrix and vector notation to get, 

9 = (Pbfs. 


(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 


This expression is a relationship between the samples of g(x ) E Vo and the exact 
scaling function coefficients. In general we need a relationship between the samples 
of f(x) E L 2 (R) and an approximate set of scaling function coefficients <r, 

/= (Pb) T 9. (89) 


This completes the matrix form of the quadrature formula. The next subsection will 
derive the matrix form of the relationship between the scaling function and its dual. 

5.2 The Scaling Function and its Dual 

Recall once more the relationship between the scaling function and its dual, 

$% ,j ^ • (90) 

In terms of the B-spline expansions we get, 

Sij = f b m {x)bi(x)dx. (91) 

m l J 
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Let B denote the matrix form of the above integral. That is, B contains the projec- 
tions of the B-splines projected onto other B-splines. The above equation becomes, 

= < 92 ) 

l m 

Now, rename the expression in the brackets 7 to get, 

ki = 5Z r j,/7.y- ( 93 ) 

/ 

Transpose r to get, 

s ij = Yhn r L- ( 94 ) 

/ 

The above equation has the following matrix form: 

/ = PBR t (95) 

This is the matrix form of the requirement that the scaling function be orthogonal 
to its dual under translation. The next subsection will derive the derivative projection 
matrix. 


5.3 The Derivative Projection Matrix 

Recall from an earlier section that the elements of the derivative projection matrix 
come from the following inner product: 

dj,i =< > • ( 96 ) 


Replace the scaling functions with their B-spline expansions to get, 

d 'j t i — J 


m t l 


Let the matrix A contain the results of the above integration, 


Now, group the terms as, 


dj y i — pt,m r jjAl )Tn . 

mj 


dj,» — E T i,lAkm\, 

m l 


(97) 


(98) 


( 99 ) 
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and rename the bracketed expression as 7 to get, 


dj,i — Pi,mHj,m • (100) 

m 

Transpose p to get, 

d i,i = (101) 

m 

and the matrix form of this equation is, 

D = RAP t . (102) 

This completes the derivation of the derivative projection matrix. The following 
subsection will combine the results from all the subsections in this section to get the 
matrix form of the differentiation matrix. 

5.4 The Differentiation Matrix 

The differentiation matrix is, 

V = CDC~\ (103) 

where 

C = ( Pb) T . (104) 

From above we found that, 

D = RAP t , (105) 

and 

I = PBR t . (106) 

From equation (106) we get, 

I = RB t P t . (107) 

Solve for R and use the fact the B is symmetric to get, 

R = (BP t )-\ (108) 
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Use this expression in equation (105) to get, 

D = (F r ) _1 B~ x AP t . (109) 

Now, using the expression for the quadrature we get the differentiation matrix, 

D = b T P T (P T )-'B- l AP T {P T )-\b T )-\ (110) 

or 

V = b T B- 1 A(b T )- 1 . (Ill) 

Let me restate two notational definitions. First, the elements of the matrix ‘6’ are 
b m ,k which denotes sampling basis ‘m’ at position i k'. Second, the elements of the 
matrix ‘ B ’ are / b m (x)bk(x)dx. 

Note that in the case of periodic boundary conditions that all three of the ma- 
trices 6, A, and B are circulant. They, therefore, commute yielding B~ l A as the 
differentiation matrix. 
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6 Examples with Periodic Boundary Conditions 

In this section explicit examples for the 1st and 3rd order spline-based wavelet differ- 
entiation matrices will be given. First, the accuracy will be illustrated explicitly by 
calculating the semi-discrete Fourier transforms for the filters a, b, and d. Then the 
matrix forms, A, B, and D will be given. 

6.1 The Differentiation Matrix for a 1-st Order B-Spline 

The first order spline is the convolution of two box functions. It is defined as fd l {x) = 
x + 1 for x 6 [—1,0], /? J (x) = 1 — x for x 6 (0, 1] and 0 otherwise. 

6.1.1 The Filters a, b , and d 

For the 1-st order B-spline the filters a and b are, 



a = {-1/2,0, 1/2}, 

(112) 

and 

^ = {1/6, 2/3, 1/6}. 

(113) 

The filter d is defined such that, 



b 3 * d = a. 

(114) 

Therefore, 

II 

(115) 

where, 

4 e* 

II 

Ms 

^3 

(116) 

and, 

&{() = 1/3(2 + cos(£)). 

(117) 

d(£) will have the form, 

d(0 = + c£ q + hot, 

(118) 
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and must satisfy 


« + cO 1/3(2 + cos(O) = •sin({). 


(119) 


Consider the first few terms of each series, 


(i£ + c£ 9 )l/3(2 + 1 - Cl 2 + Cl 24 + hot) = t(£ - Cl 6 + Cl 120 + hot), (120) 


or 


(i£ + cC)( 1 - Cl 6 + Cl™ + hot) = i{( - Cl 6 + £7120 + hot), (121) 

and finally, 

t(£ - CP) + iCp2 + c£* + hot = i(t - CP) + *'£7120 + hot, (122) 
with the conclusion that 

q = 5. (123) 

This is exactly the requirement for 4-th order accuracy. 


6.1.2 The Differentiation Matrix 


In matrix form, the filter a is, 



0 

1/2 

0 

0 

0 


0 

0 

0 


0 

1/2 


- 1/2 

0 

1/2 

0 

0 


0 

0 

0 


0 

0 


0 

- 1/2 

0 

1/2 

0 


0 

0 

0 


0 

0 


0 

0 

- 1/2 

0 

1/2 


0 

0 

0 


0 

0 

A 

0 

0 

0 

- 1/2 

0 

1/2 

0 

0 


0 

0 

A — 

0 

0 

0 

0 

- 1/2 

0 

1/2 

0 


0 

0 


0 

0 

0 

0 

0 

- 

1/2 

0 

1/2 


0 

0 


0 

0 

0 

0 

0 


0 

- 1/2 

0 

1/2 

0 


0 

0 

0 

0 

0 


0 

0 

- 1/2 


o 1/2 


. 1/2 

0 

0 

0 

0 


0 

0 

0 

— 

1/2 

o J 

and the filter b 

is, 














' 2/3 1/6 

0 

0 

0 

0 

0 

0 

0 

1/6 ' 




1/6 2/3 

1/6 

0 

0 

0 

0 

0 

0 

0 




0 

1/6 

2/3 

1/6 

0 

0 

0 

0 

0 

0 




0 

0 

1/6 

2/3 

1/6 

0 

0 

0 

0 

0 


D 


0 

0 

0 

1/6 

2/3 

1/6 

0 

0 

0 

0 




0 

0 

0 

0 

1/6 

2/3 

1/6 

0 

0 

0 




0 

0 

0 

0 

0 

1/6 

2/3 

1/6 

0 

0 




0 

0 

0 

0 

0 

0 

1/6 

2/3 

1/6 

0 




0 

0 

0 

0 

0 

0 

0 

1/6 

2/3 

1/6 




L 1/6 0 

0 

0 

0 

0 

0 

0 

1/6 

2/3 



(124) 


( 125 ) 
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The derivative projection matrix D, which recall is the same as the differentiation 


matrix T> for periodic boundary conditions, is 

D = B~ X A. (126) 

A 10 by 10 version of the derivative projection matrix for the current scenario is, 


0 

.804 

-.215 

.057 

-.014 

0 

.014 

-.057 

.215 

-.804 

-.804 

0 

.804 

-.215 

.057 

-.014 

0 

.014 

-.057 

.215 

.215 

-.804 

0 

.804 

-.215 

.057 

-.014 

0 

.014 

-.057 

-.057 

.215 

-.804 

0 

.804 

-.215 

.057 

-.014 

0 

.014 

.014 

-.057 

.215 

-.804 

0 

.804 

-.215 

.057 

-.014 

0 

0 

' .014 

-.057 

.215 

-.804 

0 

.804 

-.215 

.057 

-.014 

-.014 

0 

.014 

-.057 

.215 

-.804 

0 

.804 

-.215 

.057 

.057 

-.014 

0 

.014 

-.057 

.215 

-.804 

0 

.804 

-.215 

-.215 

.057 

-.014 

0 

.014 

-.057 

.215 

-.804 

0 

.804 

.804 

-.215 

.057 

-.014 

0 

.014 

-.057 

.215 

-.804 

0 


(127) 

6.2 The Differentiation Matrix for a 3-rd Order B-Spline 

The 3-rd order B-spline is the convolution of two 1-st order B-splines. The (non- 
normalized) 3-rd order B-spline is defined as, 0 3 (x) = (a: + 2) 3 for x € [-2,-1), 
fi 3 {x) = 4 — 6a: 2 — 3a: 3 for x € [—1,0), fi 3 {x) = 4 — 6a: 2 + 3a; 3 for x € [0,1), 
fl 3 (x) = (2 — x) 3 for x € [1,2], and 0 otherwise. 

6.2.1 The Filters a, b , and d 

The explicitly calculated filters for the 3-rd order B-spline are, 

a = {-1/20,-14/5,-49/4,0,49/4,14/5,1/20}, (128) 


b 7 = {1/140,6/7,1.191/140,604/35,1191/140,6/7,1/140}. (129) 

Normalization is not necessary since the filter a is divided by the filter b 7 . The 
semi-discrete Fourier transforms are, 

a(0 = »(49/2sin(f) + 28/5sin(2£) + l/10sin(3£)), (130) 
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and 

b 7 (£) = 604/35 + 1191/70 cos(£) + 12/7cos(2£) + l/70cos(30- (131) 

Again, the semi-discrete Fourier transforms are related by, 

i 7 U)d(() = i(0- (132) 

In a manner analogous to the previous subsection, we equate powers in the Taylor 
series of the semi-discrete Fourier transforms and come to the conclusion that, 

d(£) = i{ + c? + hot. (133) 

This is precisely the condition that is required for 8-th order accuracy. 

6.2.2 The Differentiation Matrix 

The matrices for the filters a and b are, 

A= (134) 


0 

49/4 

14/5 

1/20 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

- 49/4 

0 

49/4 

14/5 

1/20 

0 

0 

0 

- 1/20 

- 14/5 

- 14/5 

- 49/4 

0 

49/4 

14/5 

1/20 

0 

0 

0 

- 1/20 

- 1/20 

- 14/5 

- 49/4 

0 

49/4 

14/5 

1/20 

0 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

0 

49/4 

14/5 

1/20 

0 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

0 

49/4 

14/5 

1/20 

0 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

0 

49/4 

14/5 

1/20 

1/20 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

0 

49/4 

14/5 

14/5 

1/20 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

0 

49/4 

49/4 

14/5 

1/20 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

0 


and 


B = (135) 


604/35 

1191/140 

6/7 

1/140 

0 

0 

1/140 

6/7 

1191/140 

1191/140 

604/35 

1191/140 

6/7 

1/140 

0 

0 

1/140 

6/7 

6/7 

1191/140 

604/35 

1191/140 

6/7 

1/140 

0 

0 

1/140 

1/140 

6/7 

1191/140 

604/35 

1191/140 

6/7 

1/140 

0 

0 

0 

1/140 

6/7 

1191/140 

604/35 

1191/140 

6/7 

1/140 

0 

0 

0 

1/140 

6/7 

1191/140 

604/35 

1191/140 

6/7 

1/140 

1/140 

0 

0 

1/140 

6/7 

1191/140 

604/35 

1191/140 

6/7 

6/7 

1/140 

0 

0 

1/140 

6/7 

1191/140 

604/35 

1191/140 

1191/140 

6/7 

1/140 

0 

0 

1/140 

6/7 

1191/140 

604/35 
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The derivative projection matrix is 


D = B~ X A, 


D = (136) 


0 

.938 

-.387 

.184 

-.078 

0 

.078 

-.184 

.387 

-.938 

-.938 

0 

.938 

-.387 

.184 

oo 

o 

f 

0 

.078 

-.184 

.387 

.387 

-.938 

0 

.938 

-.387 

.184 

-.078 

0 

.078 

-.184 

-.184 

.387 

-.938 

0 

.938 

-.387 

.184 

-.078 

0 

.078 

.078 

-.184 

.387 

-.938 

0 

.938 

-.387 

.184 

-.078 

0 

0 

.078 

-.184 

.387 

-.938 

0 

.938 

-.387 

.184 

-.078 

-.078 

0 

.078 

—.184 

.387 

-.938 

0 

.938 

-.387 

.184 

.184 

-.078 

0 

.078 

-.184 

.387 

-.938 

0 

.938 

-.387 

-.387 

.184 

-.078 

0 

.078 

-.184 

.387 

-.938 

0 

.938 

.938 

-.387 

.184 

-.078 

0 

.078 

-.184 

.387 

-.938 

0 
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7 Nonperiodic Boundary Conditions 

In this section the differentiation matrix will be given for one particular boundary 
construction. The construction considered here is one of truncated B-splines. That 
is, away from the boundaries one simply shifts the B-splines by one to generate a 
spline basis. At the boundaries we simply continue this shifting. If part of the B- 
spline goes beyond the domain then simply truncate. This construction maintains 
the approximation properties of the spline subspace. That is, any n-th order spline 
can be generated across the domain using n-th order B-splines truncated in this way. 


7.1 1st Order Spline 


Using the boundary construction outlined in the previous paragraph and the 1st order 
B-spline one gets the following matrices: 



' -1/2 

1/2 

0 

0 

0 

0 

0 

0 

0 

0 


-1/2 

0 

1/2 

0 

0 

0 

0 

0 

0 

0 


0 

-1/2 

0 

1/2 

0 

0 

0 

0 

0 

0 


0 

0 

-1/2 

0 

1/2 

0 

0 

0 

0 

0 

A 

0 

0 

0 

-1/2 

0 

1/2 

0 

0 

0 

0 

A = 

0 

0 

0 

0 

-1/2 

0 

1/2 

0 

0 

0 


0 

0 

0 

0 

0 

-1/2 

0 

1/2 

0 

0 


0 

0 

0 

0 

0 

0 

-1/2 

0 

1/2 

0 


0 

0 

0 

0 

0 

0 

0 

-1/2 

0 

1/2 


0 

0 

0 

0 

0 

0 

0 

0 

-1/2 

1/2 


and the matrix B is, 


1/3 

1/6 

0 

0 

0 

0 

0 

0 

0 

0 

1/6 

2/3 

1/6 

0 

0 

0 

0 

0 

0 

0 

0 

1/6 

2/3 

1/6 

0 

0 

0 

0 

0 

0 

0 

0 

1/6 

2/3 

1/6 

0 

0 

0 

0 

0 

0 

0 

0 

1/6 

2/3 

1/6 

0 

0 

0 

0 

0 

0 

0 

0 

1/6 

2/3 

1/6 

0 

0 

0 

0 

0 

0 

0 

0 

1/6 

2/3 

1/6 

0 

0 

0 

0 

0 

0 

0 

0 

1/6 

2/3 

1/6 

0 

0 

0 

0 

0 

0 

0 

0 

1/6 

2/3 

1/6 

0 

0 

0 

0 

0 

0 

0 

0 

1/6 

1/3 


Previously it was shown that the differentiation matrix has the form, 


(137) 


(138) 


V = b T B~ 1 A(b T )-\ 


(139) 
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Table 1: Calculated order of accuracy for a 1st order spline- based wavelet differenti- 
ation matrix for truncated B-spline boundary construction. 


For n = 1 the matrix b is the identity. 


7.1.1 Computed Accuracy 

Again, one can estimate the order of accuracy by applying the above defined matrix 
T> to evenly-spaced samples of the function 

1 

2 + sin(x) ’ 

and comparing to the exact derivative. 

In table (7.1.1) ‘Inner Error’ denotes the error at the middle 50% of the grid 
points, and the boundary error is the error at the two outermost grid points at each 
end of the interval. From the table it can be seen that the superconvergence is lost at 
the boundaries. That is, in the middle of the interval the 4th-order superconvergence 
is obtained, whereas at the boundary the differentiation is only lst-order accurate. 

7.2 Example for 3-rd Order Spline 

For the 3rd order B-splines the calculations are similarly straightforward. There is, 
however, one problem: the matrix b is ill-conditioned. We need the inverse of b and 
for grids larger than about 32 points the errors introduced by inverting b make it 
very difficult to determine the accuracy. One can avoid this problem when testing 
accuracy by simply not including b in the calculations. That is, b comes from the 
quadrature equations and plays no role in the differencing. The following are the 
matrices A and B produced by the truncated B-spline boundary construction: 
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A= (141) 


- 1/2 

- 9/20 

9/10 

1/20 

0 

0 

0 

0 

0 

0 

- 71/20 

-8 

87/10 

14/5 

1/20 

0 

0 

0 

0 

0 

- 19/10 

- 127/10 

- 1/2 

49/4 

14/5 

1/20 

0 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

0 

49/4 

14/5 

1/20 

0 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

0 

49/4 

14/5 

1/20 

0 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

0 

49/4 

14/5 

1/20 

0 

0 

0 

0 

- 1/20 

- 14/5 

- 49/4 

0 

49/4 

14/5 

1/20 

0 

0 

0 

0 

- 1/20 

- 14/5 ■ 

- 49/4 

1/2 

127/10 

19/10 

0 

0 

0 

0 

0 

- 1/20 ■ 

- 14/5 

- 87/10 8 

71/20 

0 

0 

0 

0 

0 

0 

- 1/20 

- 9/10 9/20 

1/2 





B = 





( 

1/7 

129/140 

3/7 

1/140 

0 

0 


0 

0 

0 

129/140 

302/35 

531/70 

6/7 

1/140 

0 


0 

0 

0 

3/7 

531/70 

599/35 

1191/140 

6/7 

1/140 

0 

0 

0 

1/140 

6/7 

1191/140 

604/35 

1191/140 6/7 

1/140 

0 

0 

0 

1/140 

6/7 

1191/140 

604/35 

1191/140 

6/7 

1/140 

0 

0 

0 

1/140 

6/7 

1191/140 604/35 1191/140 

6/7 1/140 

0 

0 

0 

1/140 

6/7 

1191/140 599/35 

531/70 

3/7 

0 

0 

0 

0 

1/140 

6/7 

531/70 

302/35 129/140 

0 

0 

0 

0 

0 

1/140 

3/7 

129/140 

1/7 J 


7.2.1 Computed Accuracy 


If we apply the matrix T> = B 1 A to the evenly-spaced samples of the function 

1 

2 + sin(£r) ’ 

for a fixed wave number £ it is difficult to determine the accuracy both at the bound- 
aries and in the middle of the interval at the same time. That is, it was proven that 
the accuracy ‘in the middle’ is 2n + 2 or 8 when n = 3. For the current boundary 
construction, however, the accuracy is 3 at the boundaries. This vast difference in 
accuracy requires that the frequency of the test function be low in order to test the 
boundary accuracy and high in order to test the superconvergence in the middle of 
the interval. 

In table (7.2.1) the wave number is £ = 27 r. This wave number is sufficiently high 
to test the superconvergence. 



V = B~'A 


(144) 


29 


III] 

1 


Error 

Ratio 

Boundary 
/i Error 

Error 

Ratio 

Total 
l\ Error 

Error 

Ratio 

16 

mmm 


.022 


.0037 


32 


1944=2 10 - 9 

6.6e -4 


4.9e“ 5 


64 

2.4e~ 10 

750 = 2 95 

3.5e" 5 


— 


128 

eeemm 

324 = 2 8-3 

2.8e -6 

i^H 



256 

4.3e -15 

mojfifl 

3.1e -7 


2.7e -9 

U3BHI 


Table 2: Calculated order of accuracy for a 3rd order spline-based wavelet differen- 
tiation matrix for the truncated B-spline boundary construction. The test function 
has a high frequency in order to capture the 8th order differentiation accuracy. 


Grid 

Size 

— 

Error 

Ratio 

Boundary 
/i Error 

Error 

Ratio 

Total 
/] Error 

Error 

Ratio 

16 



1.8e -6 


WMSM 


32 

QQII 

390=2 8 ' 6 

KSSfli 

5.1=2 24 


11.7 = 2 36 

64 

8.9e -15 

4607 = 2 12 ' 2 

5.1e -8 

6.9 = 2 28 

UQSfl 

QVjBlil 

128 



6.8e -9 




256 

8.2e -15 



■fcgBIWSI 

7.4e- 12 


512 



l.le -10 

8 = 2 3 

QS| | 



Table 3: Calculated order of accuracy for a 3rd order spline-based wavelet differen- 
tiation matrix for the truncated B-spline boundary construction. The test function 
has a low frequency in order to capture the 3rd order differentiation accuracy. 

Now the lower frequency £ = 1 is used to test the accuracy at the boundaries. See 
table (7.2.1). Note that in table (7.2.1) that a few of the entries in the ‘Error Ratio’ 
columns are blank. This is done whenever the ‘Inner Error’ is of the order 10 -15 and 
can no longer be resolved on the computer. Table (7.2.1) indicates the accuracy only 
at the boundaries. As in table (7.2.1) it can be seen that the boundary accuracy is 
3rd order and, consequently, the superconvergence is lost at the boundaries. 

The next section contains a few comments on this loss of superconvergence. 
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8 Loss of Super convergence 


In [5] it was shown that under a given boundary construction the Daubechies-based 
wavelet superconvergence which was proved in [4] is lost. In this paper we have 
encountered a similar phenomenon. I believe that this in not simply a function of 
the the particular boundary functions which were chosen but a general characteristic 
for numerical methods which are based on approximation by translation. In fact, 
Gottlieb, et. al, see [3], have shown that the superconvergence encountered in the 
finite element method under periodic boundary conditions cannot be recovered when 
nonperiodic boundary conditions are present: ‘... no matter how the approximating 
subspace S h is chosen, the superconvergence property is lost if there are characteristics 
leaving the domain.’ 

Superconvergence is a property which depends on the overlap of bases functions 
with neighboring bases functions. At the boundaries the overlap only occurs on one 
side. For a Daubechies-based wavelet system where the bases functions are orthogonal 
there is a minimum amount of overlap which must occur at the boundary in order 
to maintain the superconvergence. That is, for a Daubechies 4 wavelet, 4-th order 
differentiation accuracy is encountered. In order to maintain the superconvergence 
one must have at least 3rd order differentiation accuracy at the boundaries. That is, 
four conditions must be satisfied requiring an overlap of four basis functions at the 
boundary. As of now, there are no boundary constructions satisfying this criterion. 

For spline-based wavelet systems the situation is similar but slightly more compli- 
cated since the underlying functions, B-splines, are not orthogonal. Again, however, 
a general rule-of-thumb is that the overlap of the B-splines at the boundary should 
be the same as the superconvergence encountered with periodic boundary conditions. 
That is, for the 1st order B-spline we have shown here that the differentiation is 4th 
order accurate where periodic boundary conditions are assumed. This would require 
an overlap of 4 boundary functions in order to have 3rd order boundary accuracy 
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which maintains the 4th order overall accuracy. 
The next section concludes this paper. 



9 Conclusion 


Let me restate a few of the salient features of this paper. First of all, circularity plays 
a prominent role when periodic boundary conditions are assumed. The benefit of 
working with circular matrices is that they commute. In the periodic case this implies 
that the differentiation matrix is the same as the derivative projection matrix, 

V = C~'DC = D. 


The proof of superconvergence is, however, the most meaningful result of the pa- 
per. I have shown that the matrix D displays accuracy of order 2n + 2 where n is 
the order of the underlying B-spline used to construct the scaling functions. Also, it 
was seen that this accuracy holds regardless of of whether the scaling functions are 
chosen to be orthogonal under translation or not. That is, the scaling functions can 
be required to be B-splines, C-splines, or orthogonal. The parameters which dictate 
the properties of the scaling function ‘divide out’ when the differentiation matrix is 
calculated. Because of this feature, only one theory is necessary to illustrate the ac- 
curacy for the spline-based differentiation matrix when periodic boundary conditions 
are assumed. 

When the boundary conditions are, on the other hand, not periodic it has been 
shown that for boundary functions constructed from truncated B-splines that the su- 
perconvergence is lost at the boundary. I speculate that this loss of superconvergence 
will continue to be a characteristic of wavelet bases defined on an interval. 

It is, also, interesting to note the relationship between spline-based wavelet meth- 
ods and compact finite difference methods. That is, the underlying character of a 
spline-based wavelet method can be understood in terms of compact finite difference 
schemes. 


33 



References 


[1] I. Daubechies, “Orthonormal Basis of Compactly Supported Wavelets”, Comm. 
Pure Appl. Math., 41 (1988), pp. 909-996. 

[2] P. Davis, “Circulant Matrices”, Wiley-Interscience, 1979. 

[3] David Gottlieb, Bertil Gustafsson, Pelle Olsson, and Bo Strand, “On the Super- 
convergence of Galerkin Methods for Hyperbolic IB VP”, July 29, 1993, submit- 
ted to SIAM J. Sci. Comp. 

[4] L. Jameson, “On The Wavelet Based Differentiation Matrix”, Journal of Scien- 
tific Computing, September 1993. 

[5] L. Jameson, “Daubechies-Based Wavelets on an Interval”, to be submitted to 
SIAM J. Sci. Comp. 

[6] J. Liandrat and P. Tchamitchian, “Resolution of the ID Regularized Burgers 
Equation using a Spatial Wavelet Approximation Algorithm and Numerical Re- 
sults”, ICASE Report No. 90-83 (1990). 

[7] I. J. Schoenberg, “Contribution to the problem of approximation of equidistant 
data by analytic functions”, Quart. Appl. Math. 4 (1946), 45-99, 112-141. 

[8] I. J. Schoenberg, “Cardinal Spline Interpolation”, CBMS-NSF series in Applied 
Math 12, SIAM Publ., Philadelphia, PA, 1973. 

[9] G. Strang, “Wavelets and Dilation Equations: A Brief Introduction”, SIAM 
Review, Vol. 31, No. 4, pp. 614-627, Dec. 1989. 

[10] M. Unser and A. Aldroubi, “Polynomial Splines and Wavelets- A Signal Process- 
ing Perspective”, Wavelets- A Tutorial in Theory and Applications, C. K. Chui 
(ed.) pp. 91-112 (1992), Academic Press, Inc. 


34 



[11] M. Unser, A. Aldroubi, and M. Eden, “The L 2 polynomial spline pyramid: a dis- 
crete multiresolution representation of continuous signals”, IEEE Pattern Anal. 
Mach. Intell., to appear. 

[12] Wei Cai and JianZhong Wang, “Adaptive Wavelet Collocation Methods for Ini- 
tial Value Boundary Problems of Nonlinear PDE’s”, ICASE Report No. 93-48 
NASA CR-191507, and submitted to SIAM J.Num. Anal. 


35 




REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect ol this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports 1215 Jefferson 
Davis Highway Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington. DC 20503 

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 

November 1993 

3. REPORT TYPE AND DATES COVERED 

Contractor Report 

«. TITLE AND SUBTITLE 

ON THE SPLINE-BASED WAVELET 
DIFFERENTIATION MATRIX 

5. FUNDING NUMBERS 

C NAS1-19480 
WU 505-90-52-01 

6. AUTHOR(S) 

Leland Jameson 


7. PERFORMING ORGANIZATION NAME(S) AND AODRESS(ES) 

Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23681-0001 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

ICASE Report No. 93-80 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23681-0001 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA CR-191557 
ICASE Report No. 93-80 

11. SUPPLEMENTARY NOTES 

Langley Technical Monitor: Michael F. Card 
Final Report 

To be submitted to SIAM Journal on Numerical Analysis 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 

U nclassified-U nlimited 




Subject Category 64 | 

13. ABSTRACT (Maximum 200 words) , 

The differentiation matrix for a spline-based wavelet basis will be constructed. Given an n-th order spline basis it 
will be proven that the differentiation matrix is accurate of order 2n + 2 when periodic boundary conditions are 
assumed. This high accuracy, or superconvergence, is lost when the boundary conditions are no longer periodic. 
Furthermore, it will be shown that spline-based bases generate a class of compact finite difference schemes. 


16. PRICE CODE 

A03 

17. SECURITY CLASSIFICATION 10. SECURITY CLASSIFICATION! 19. SECURITY CLASSIFICATION] 20. LIMITATION 

OF REPORT OF THIS PAGE OF ABSTRACT I OF ABSTRACT 

Unclassified Unclassified I 


14. SUBJECT TERMS 

differentiation matrix; wavelets; superconvergence 


15. NUMBER OF PAGES 

37 


NSN 7540-01-280-5501 


☆ U.S. GOVERNMENT PRINTING OFFICE: 1993 - S2H-064/S60M 


Standard Form 298(Rev. 2-89 

Prescribed by ANSI Std 239-18 
298-102 













